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The main point raised in the Comment of Huopaniemi 
et al. [H] concerns the scaling of the mean time (r u ) it 
takes a polymer of length N, threaded halfway in a nar- 
row pore, to unthread, in the absence of any external field 
or pulling force on the polymer (i.e., for unbiased trans- 
location). As argued in our paper the mean dwell 
time (Td) that a translocating polymer spends in the pore 
scales with polymer length TV in the same way as the un- 
threading time. We specifically studied this for polymers 
in three dimensions, whose dynamics is described by the 
combination of reptation and Rouse dynamics; i.e., hy- 
drodynamics is neglected. 

On the theoretical side, a relevant time scale for this 
problem is the Rouse time tr, which is the longest time 
scale for a polymer in bulk solution to relax in the ab- 
sence of external forces. As a function of polymer length 
N, the Rouse time scales as tr ~ N 1+2u . We verified 
that this scaling holds in the lattice polymer used in our 
simulations, and also that the Rouse time is the longest 
time scale for a polymer tethered to a fixed membrane 
0. Since the mobility of a polymer threaded in a pore 
will not exceed that of an unrestricted polymer, it follows 
that (Td) > tr @, In the existing literature, there is 
no theoretical argument for this inequality to reduce to 
an equality. There is however numerical evidence in 2D 
that is inequality is saturated in Ref. jsj, as well as in 
the Comment and the earlier works of the authors of the 
Comment. With the number s of the monomer located in 
the pore taken as a reaction coordinate, a consequence of 
the above inequality is that the diffusion of this reaction 
coordinate has to be anomalous, i.e. the mean squared 
displacement (As 2 (i)) ~ t a with an anomalous dynamics 
exponent a < 2/(1 + 2v). Again, there is no theoretical 
argument why also this inequality should be saturated. 

Since both a and the scaling (t c i) ~ tr in Ref. [2] were 
obtained solely from a single set of simulations to calcu- 
late (Td), our first remark concerns the factual misrep- 
resentations by the Comment's authors, to suggest that 
the results of Ref. [2] are "well-established". The authors 
of the comment wished to settle these with simulations 
alone. 

There is plenty of numerical evidence that points to- 
wards different scaling of (Td) than tr, both in 3D and 
in 2D. In our paper we reported a numerical exponent 
2.40 ±0.05 for unbiased translocation in 3D [5[. Another 
group, using a completely different polymer model, re- 
ported an exponent 2.52 ±0.04 Q. In subsequent works, 



we have provided a full theoretical description of this 
problem, leading to the result (to) ~ N 2+v both in 3D 
[fl@] and in 2D [§]. This theoretical description is sup- 
ported by high-precision numerical simulations, for which 
the 3D results we provide below in Table I. 
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100 


65136 


0.434 


150 


183423 


0.428 


200 


393245 


0.436 


250 


714619 


0.445 


300 


1133948 


0.440 


400 


2369379 


0.437 


500 


4160669 


0.431 



Table I: Median unthreading time over 1,024 runs for each value 
of the polymer length N in 3D. Data taken from Ref. Q|. 

The only numerical evidence contradicting our theory , 
as far as we are aware of, is the newly produced numerical 
result in the Comment, and that of Wei et al. [10]. How- 
ever, it is unclear whether the simulations using GRO- 
MACS by the Comment's authors or those used by Wei et 
al. implement purely Rouse and reptation dynamics, as 
is the assumption in all theoretical work mentioned here. 
Moreover, the authors [3, [T3| have taken N only up to 
200, from which an attempt to recover scaling results for 
(Td), in our opinion, is misleading. Specially, since one 
is dealing with a numerical difference of order 10%, the 
discrepancy between different simulation results can eas- 
ily be due to finite- N effects; replacing N by N + y/N or 
N — \/N produces double-logarithmic plots in which the 
data can be fitted about equally well by straight lines, 
however with exponents that deviate easily 10% or more. 
We do not believe that this apparent discrepancy can be 
resolved by simulations alone. 

Although the full derivation of the result (Td) ~ N 2+y 
can be found elsewhere @, @, [§] , for the sake of complete- 
ness we summarize it below. 

Translocation takes place via the exchange of mono- 
mers through the pore. This exchange responds to (j)(t), 
the difference in chain tension perpendicular to the mem- 
brane; simultaneously, 4>(t) adjusts to v(t) = s(t), the 
transport velocity of monomers across the pore, as well! 
With As(£) as the total number of monomers translo- 
cated from one side to the other in the time interval [0, t], 
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and 4>{t) playing the role of chemical potential difference 
across the pore, the two variables As(i) and 4>(t) are 
conjugate to each other in the thermodynamic sense. In 
the presence of memory effects, they are related to each 
other by (j>{t) = J dt' fj.(t-t')v(t') via the memory kernel 
fJ,(t), which can be thought of as the (time-dependent) 
'impedance' of the system. This relation can be inverted 
to obtain v(t) = j Q dt'a(t — t')(p(t'), where a(t) can be 
thought of as the 'admittance'. In other words, in the 
Laplace transform language, fi(k) — a^ 1 (k), where k 
is the Laplace variable representing inverse time. Ad- 
ditionally, via the fluctuation-dissipation theorem, they 
are related to the respective autocorrelation functions as 
H(t - t') = (0(i)0(f ))»=o and a(t - t') = («(t)«(f )}#=o- 

In Ref. f§] we showed that fi{t) ~ t~~^ exp(-t/r R ) 
in 3D. This implies that the translocation dynamics is 
anomalous for t < tr, i.e., (As 2 (t)) = J Q dt'(t — t')a(t'), 
the mean-square displacement of the monomers through 
the pore, behaves as t ai for some a% < 1. Beyond the 
Rouse time the translocation dynamics becomes simply 
diffusive. From the behaviour of above, it is easily 
shown that ot\ = l+%v ' having ignored the exp(— i/r^) 
term for t < tr, one obtains ji(k) ~ k~ x + 2 » , implying 
a(k) ~ fciw, i.e., a(t) ~ rw ] which yields cti = jtt^. 

Thus, for t < tr, (As 2 {t)) - and for f > t r 



(As 2 {t)) ~ t, which together yield (r d ) - iV 2+l/ , both 
in 3D and 2D. Moreover, using high-precision simulation 
data, we demonstrated that in 2D the probability dis- 
tribution of the dwell time P{Td), behaves as P(rd) ~ 
V{T d /N 2+v )/N 2+v , with a scaling function V(t) [|. 

To conclude, to date no theoretical argument has been 
reported for why (r^) should scale as tr. In fact there 
is a derivation why the scaling with polymer length N 
for the two should differ @, U 0| • Numerical evidence 
0, 0, B [§] also points towards a scaling of (r^) different 
from Tfl, apart from those due to the Comment's authors 
and that due to Wei et al. [Ifj. The theoretical formalism 
that yields (r^) ~ N 2+l> for unbiased translocation also 
works beautifully for translocation mediated by a pulling 
force at the head of the polymer and field-driven 
translocation [jl. fl2l|. providing a solid unified theoretical 
understanding of the dynamics of translocation, based on 
the well-known laws of polymer physics. 

We end our reply with the additional observation that 
the expression for (r^) for pore length L provided in the 
Comment [J, Hj> too, is incorrect. There is general agree- 
ment that the monomers inside the pore show anoma- 
lous diffusion with some exponent a < 2/(1 + 2u), as 
discussed above. With this kind of dynamics, the time 
to travel over a distance L has to increase faster than 
quadratically. 
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